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We present some new derivations of the effect of a plane gravitational wave on a light ray. A simple 
interpretation of the results is that a gravitational wave causes a phase modulation of electromagnetic 
waves. We arrive at this picture from two contrasting directions, namely null geodesics and Maxwell’s 
equations, or geometric and wave optics. Under geometric optics, we express the geodesic equations 
in Hamiltonian form and solve perturbatively for the effect of gravitational waves. We find that the 
well-known time-delay formula for light generalizes trivially to massive particles. We also recover, 
by way of a Hamilton-Jacobi equation, the phase modulation obtained under wave optics. Turning 
then to wave optics - rather than solving Maxwell’s equations directly for the fields, as in most 
previous approaches - we derive a perturbed wave equation (perturbed by the gravitational wave) 
for the electromagnetic four-potential. From this wave equation it follows that the four-potential 
and the electric and magnetic fields all experience the same phase modulation. Applying such a 
phase modulation to a superposition of plane waves corresponding to a Gaussian wave packet leads 
to time delays. 


I. INTRODUCTION 

While the existence of gravitational waves is secure 
from their indirect consequences, namely, the slow orbital 
infall of binary pulsars BE], their direct observation is 
nonetheless an eagerly awaited event. 

Gravitational-wave detectors seek to measure the effect 
of light traveling through a “waving” space-time. Inter¬ 
ferometers with arms that are kilometers long have been 
in operation for some time [3] and, with recent increases 
in sensitivities mu it seems likely that detection of waves 
originating from binary coalescence of stellar-mass black 
holes or neutron stars may be around the corner. Such 
interferometers are commonly introduced mu through a 
picture of emitter and receiver as freely falling test parti¬ 
cles experiencing a time-dependent acceleration from the 
gravitational wave, with an interferometric setup measur¬ 
ing the resulting changes in optical paths. This simple 
picture is applicable if the light path is much shorter 
than the wavelength of the gravitational wave, a rea¬ 
sonable approximation for ground-based interferometers. 
The situation may be very different for a space-based in¬ 
terferometer mm, or an experiment with spacecraft 
clocks mm, because these could operate in the free 
spectral range, where optical path lengths are compara¬ 
ble to the gravitational wavelength. And for a pulsar tim¬ 
ing array, consisting of a selection of millisecond pulsars 
under near-continuous observation, the pulses would have 
traveled through many gravitational wavelengths before 
reaching the observer [18] . 

For a more general picture it is necessary to con¬ 
sider light propagation itself through a gravitational wave 
space-time. There are two contrasting strategies for do¬ 
ing this, which may be described as geometrical optics 
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and wave optics. In geometrical optics one considers null 
geodesics in a waving space-time, and hence calculates a 
light travel time. The basic time-delay formula was first 
derived in the 1970s mm and has been revisited more 
recently [20}l22] . Some other early derivations are now 
deprecated [23] because they involve invalid shortcuts. 
The transverse-traceless gauge is standard, but a gauge- 
invariant approach, which expresses the time delay as an 
integral of the Riemann curvature along the photon path, 
is also possible [24]. 

The use of geometrical optics, even where one wants 
to measure phase changes, is commonly justified as the 
eikonal approximation to electromagnetic waves m • 
Nonetheless, it is also possible to solve Maxwell’s equa¬ 
tions in a waving space-time [25427] . 

This work attempts to synthesize the geometrical- 
optics and wave-optics views, and to provide an intu¬ 
itive description of what a gravitational wave does to an 
electromagnetic wave. Such a comparison has previously 
been done for the small-antenna regime, that is, when 
the gravitational-wave phase does not change while light 
is traveling through the interferometer [28 . In space- 
based interferometers, however, the antenna length (say 
L ) may be comparable to the gravitational wavelength 
2ir/k , and in pulsar timing kL 1. It is thus interesting 
to consider the free spectral range. 

In Sec. m we consider photons moving on null 
geodesics. A four-dimensional Hamiltonian formulation, 
with the affine parameter being the independent vari¬ 
able and time being a coordinate, is well suited to this 
problem. We derive geodesics using two different pertur¬ 
bative approaches. In addition to the well-known time- 
delay formula, two interesting results emerge from these 
calculations: (i) a generalization to massive particles and 
(ii) the notion of a phase even in geometrical optics. 

In Sec.[lV]we consider the electromagnetic wave equa¬ 
tion through a Minkowskian space-time perturbed by a 
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FIG. 2. Following from Eq. |l]), the red curve shows the 
maximum signal possible over all inclinations, #. One does 
not gain by increasing the detector length beyond the gravi¬ 
tational wave’s wavelength. The blue curve shows the incli¬ 
nation at which this maximum occurs. 


FIG. 1. This polar plot shows the angular dependence, #, be¬ 
tween the gravitational and electromagnetic wave directions 
- the normalized form of equation ([Tj) . The radial direction is 
normalized. When the detector is short compared to the grav¬ 
itational wave period, the phase shift’s directional dependence 
reduces to sin 2 9 . When the detector is long (or the gravita¬ 
tional wave short), a more parallel alignment is favored, and 
an antiparallel disfavored. The maxima of the amplitudes are 
shown in Fig. [2] 


gravitational plane wave, and we write down a simple 
perturbed wave equation for the four-potential. The per¬ 
turbation terms, coupling the gravitational wave to the 
electromagnetic field, cause a modulation of the phase 
velocity along the light path. 

We show that the net effect of the gravitational wave 
on light is remarkably simple: the phase of the electro¬ 
magnetic (EM) wave gets modulated like 

<f> = e^- cos2 (j) (1 + cos#) sin[fc(t — r cos#)], (1) 

2k 

where we have chosen geometric units c = 1, and # and 
< j) give the propagation and polarization direction of the 
gravitational wave (GW) with respect to the EM wave 
(EMW). The phase modulation can be thought of as re¬ 
distributing the energy of the EM wave, but no net work 
is done at order e. A light pulse sent from an emitter 
to an observer can be expressed as a superposition of 
plane waves. Phase modulation alters the pulse shape, 
but the main effect is to move the pulse, thus chang¬ 
ing the arrival time. For simplicity, the above expression 
considers a gravitational plane wave with a single wave¬ 
length, one polarization state, and a particular phase. 
A real detector has to contend with a superposition of 
such contributions. Later in this paper, we will relate 


Eq. 0 to expressions in the literature on space-based 
interferometers and pulsar timing. 


II. THE METRIC 


As in most work on gravitational waves, we will work 
in the transverse-traceless gauge, with the mostly pluses 
signature convention, in which the metric 
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cos [k (z — t)] (2) 


is a weak-field solution to the vacuum gravitational field 
equations, with a gravitational wave with amplitude e 
and wave vector k traveling through it. For convenience, 
we have chosen the gravitational wave to approach from 
the —z direction, crest at the z = 0 plane, and bear only 
®-type polarization. With these choices we retain full 
generality because at any point in the calculation we are 
free to apply a spatial rotation or time shift. 

For the contravariant components of the metric, one 
needs to note that numerically h^ v = —h^ u . For future 
reference, we note the determinant of g^ u 

_ f 2 

\f^g = 1 — — cos 2 ( kz — kt) = 1 + O (e 2 ) . (3) 

We further note that the nonzero Christoffel symbols, up 
to O (e 1 ), are all the same: 

ke 

r* tx = T y ty = T* xx = Y\ y = - sin (kz - kt). (4) 
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III. GEOMETRICAL OPTICS 

In this section we consider light to be a freely falling 
trajectory, i.e., one which parallel transports its tangent 
vector. The gravitational wave will affect a photon’s mo¬ 
mentum as it travels, and thereby alter the time at which 
the receiver records its arrival. A convenient formulation 
is in terms of a Hamiltonian system [29 . The relevant 
Hamiltonian is 


H = \p^Pv9^ u (5) 

or 

H = | (~Pt + PiP 0 + p (pi ~ Pi) cos {kz - kt) (6) 

when we substitute the metric ©• The independent vari¬ 
able is the affine parameter, say A, and t is simply a co¬ 
ordinate, and p^ is the four-momentum. If H = 0 the 
trajectory corresponds to photons, whereas H < 0 de¬ 
scribes timelike geodesics. 


A. Null and timelike geodesics 


Since we are looking for the time delay, the Hamilton 
equations we need to solve are 


. d H 

t = w— = -Pt 
dpt 

• dH kc / o o \ /1 7 \ 

Pt = = -y [Py-Px) sin {kz - kt ). 


( 7 ) 


Let (T, A, Y, Z , P tl P x ,P yi Pz) be a solution with no grav¬ 
itational wave. We have 


P t = const, T = —A P t (8) 

and so on. For convenience, we set P t = —1, in effect 
choosing the scale of A. Substituting the unperturbed 
solution in the 0(e) terms in (|7|) and integrating, we ob¬ 
tain 


does not actually require the lightlike, null stipulation 
P^P* 1 = 0. In fact, the delay is valid for null and massive 
trajectories alike. Hence, the time delay © applies also 
to massive particles. For the latter, we have 


z = r cos 0 

p z = -vp t cos 0 (12) 

Px~P y = y2 Pt s i n2 0 s i n 20 

where v is the speed. The time delay 0 can be written 
as 


e v 2 sin 2 0 
2k 1 — v cos 0 


sin (kr cos 0 — kt ). 


(13) 


Taking the limit v 1 and multiplying by the uj of 
the light gives the phase change For nonrelativistic 
particles, the effect of gravitational waves on trajectories 
would be reduced by 0(v 2 ). It is tempting to contem¬ 
plate experiments with neutrinos shot through the Earth, 
with their exit times measured; but we cannot think of 
any plausible scenarios for detecting gravitational waves 
using timelike geodesics. 

Returning now to light paths, we can also rewrite the 
time delay in a form similar to well-known expressions in 
recent literature. Let r l be the unit vector along which 
a light ray travels, and let k l be the corresponding unit 
vector for the gravitational wave. Making a comparison 
with the metric < 0 , we see that the time delay © is 
nothing but 


At = 0 T -k l r l I kij ^ < ' 14 ' ) 

with the particular choice of a gravitational wave prop¬ 
agating in the positive-^ direction in one polarization 
and initial phase. The expression may be compared with 
Eq. (11) of Ref. [10! for space interferometers, and with 
Eqs. (6) and (7) of Ref. JSj for pulsar timing. 


B. A Hamilton-Jacobi formulation 


Pt 


p2 

V 


2 P x 


cos (kZ — kT) — 1. 


( 9 ) 


An alternative perturbation method for Hamiltonian 
systems is to seek a transformation from 


Integrating once again, we arrive at 


((A HfT#ik sln( ‘ z - tz)+A - 


( 10 ) 


The time delay A t = t — T = t — A, which we may express 
as 


At 


e Pi ~ Pi 
2k (p z +Pt)Pt 


sin (kz — kt) 


(ii) 


with the understanding that the unperturbed values can 
be used on the right. 

Although at the beginning of this section we referred to 
this trajectory as a model for light, the above derivation 


(x^iPn) =(t,x k ,Pt,Pk) (15) 

to a new set of canonical variables 

(XrP^) = (T.X\P t ,P k ) (16) 

such that the transformed Hamiltonian (say K) becomes 
trivially integrable. Thus 

H(p^, x») = K(P tu ) = -P? + P k P k . (17) 

We allow ourselves the expression PkPk since the trans¬ 
formed space is Euclidean. The geodesics in that space 
are 


T = -A P t X k = XP k • 


(18) 
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As before, we take the origin as the start of the geodesics. 
The canonical transformation can be specified from a 
generating function S(P /1: x fJ '), such that 


p„ = p,p 

= x^ 


dS 

t)xJ l 

dS 

OP, 


(19) 


which is not one of the textbook forms, but is easy to 
derive, by considering 


p v dx v — H dX = P u dX v -KdX 

-dS(P^x^) + d(P v x v ) - d(P v X v ) .(20) 
The relation between the Hamiltonians is then 

<2i) 

This is a form of the Hamilton-Jacobi equation, and it 
has the perturbative solution 

k ^ + ^ csh{p ^ ) - (22) 

Substituting in the Hamiltonian ([6]), it is easy to solve 
for 

s =h v tw, sia(kz ~ kt) - (23) 


Since photons have p M = huo u , this S/h equals the wave- 
optics phase modulation (40), which is the eikonal ap¬ 
proximation. 


To calculate the time delay, we take the geodesics (18) 
and apply the canonical change of variable (19). Then 


t~x* + g. 

opt 

k dS 

c k = XP k + — 
OPk 


(24) 


are geodesics in the original variables. Again, we allow 
ourselves to mix up and down indices here, because 
is Minkowskian. Now using the null condition 

P 2 t = P k P k (25) 


we have 



dS 

dpk 


(26) 


and with r = \x k \ we can simplify the above expression 
to 


-A P t = r - 


x k dS 


r dpk 


(27) 


to first order. We can now write down the difference in 
light travel time compared to the unperturbed case 


At = t — r 


dS_ 

dpt 


x k dS 
r dp k 


(28) 


which reduces to the result derived from the direct 
method. 


IV. WAVE OPTICS 

An alternative to treating light as a geodesic path is to 
consider it as an electromagnetic field excitation. We will 
proceed to write down Maxwell’s equations in the metric 
© and then find a solution for a plane electromagnetic 
wave perturbed by the gravitational wave. 


A. Maxwell’s equations 

In freely falling coordinates, Maxwell’s equations are 

X U F^ = 4t rJ^ (29) 

where the field tensor is related to a four-potential by 

F^ = d»A v - d u A‘ u . (30) 

It is understood that fields and potentials are functions of 
x To obtain the general, global form of the equations, 
the partial derivatives are promoted to covariant. This 
results in 


X U X^A U - X U X U A^ = 4t rPb (31) 

We can commute the covariant derivatives through the 
introduction of the Riemann tensor: for an arbitrary vec¬ 
tor field Vp, 


[V„,V 7 ]^-K a J?Vr. (32) 


Using this to commute the derivatives in the first term 
of (31) gives 


X^X U A U + R» a A a -X U X V A» = . (33) 

o © 


The labeled terms vanish, for the following reasons. 

O Adopting the Lorentz gauge X V A V — 0 makes this 
term fall away. 

© We make the approximation that our detector’s 
cavity is a vacuum. There are no charges in a vac¬ 
uum, so = 0; 

© We also assume there are no sources of stress, en¬ 
ergy, or momentum, and hence the Ricci tensor 
components R^ = 0. The electromagnetic wave 
itself contains some energy-momentum, of course, 
and the resulting metric perturbations can be ex¬ 
pected to be of the same order as the laser power 
in dimensionless terms. The latter quantity is 
G/c 5 x [laser power], and will be far smaller than 
even the already minuscule gravitational waves of 
interest. In other words, electromagnetic excita¬ 
tions do not induce significant gravitational ones. 
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The electromagnetic field equations reduce to 


□g^ = 0, 


(34) 


where D g = is the covariant d’Alembertian oper¬ 

ator, a global generalization of d u d u , and can be written 


= (35) 


Equation (34) are the equations of motion for the elec¬ 


tromagnetic potentials in an arbitrary gravitational back¬ 
ground, in the vacuum approximation. The next step is 
to let the gravitational wave provide the background. 


B. Electromagnetic wave equation in a 
gravitational-wave space-time 


Recalling that g = 1 to leading order O (e 1 ) [Eq. 
£ 1 ) 1 , the d’Alembertian |35l) reduces to 


= d „ [{rf v - h H d v A x ] 


rf v d v d lt -d lt h' u 'd v -h' M 'd ll d v 


=0 




(36) 


The first term is the flat space-time wave equation. 
The second term is actually zero, because when /x, z/ = 
x, y, the derivative is zero, and when /x, v — t,z, the 
metric perturbation h^ m 0. This leaves us with the 
third term as the only gravitational perturbation to the 
flat space-time result. The wave equation to solve, then, 


is 


□ 77 + 6 cos (kz — kt ) 




A" = 0. 


(37) 


The first-term box is the Minkowskian d’Alembertian. 


C. Electromagnetic plane waves 

In flat space-time, an electromagnetic plane wave can 
be expressed as 


A x = Aq exp (iuj^x u r ]^ u ), (38) 

where the wave vector uj x satisfies the null condition 
= 0. The complex constant vector Ag encodes 
the amplitude and polarization of the wave. 

It is straightforward to verify that the perturbed plane 
wave, that is, the real part of 

A x = Aq exp {iuj pii x v rf ,v + x<h) (39) 


where 


4 > = 


e - u l ■ n i *\ 

- - sm [kz — kt) 


2 k uj t 


■ UJ 


(40) 


satisfies the wave equation ( |37| ) to leading order O (e 1 ). 
This is an electromagnetic plane wave with a phase mod¬ 
ulation 4> coming from the gravitational wave. It is also 
equivalent to the geometrical optics result ( 11 ). 

To compute the components of the electric and mag¬ 
netic fields, we must evaluate the covariant derivatives. 
For the electric field, we have 


E i = pti = \yt A i _ yi A t j 

= [d l A 1 - <9* A* + + h**) T\ X A X - + h^) T\ X A X ] (41) 

= - a 1 a* - (r tx + t\ x ) a x - (r\ y + ry a«] . 


In the last line, the and h 1X1 terms have been dropped, 
because they are ~ e 2 3 . The nonzero Christoffel symbols 
are given in Similarly complicated expressions give 
the magnetic field. 

Fortunately, there is no need to work out the covari¬ 
ant derivatives in full. One just needs to consider the 
different combinations of a;, e, k that would appear in the 
covariant derivatives. At zero order, there will of course 
be terms oc uj as in flat-space plane waves. Further terms 
will be like: 

1. O(euo 2 /k ), from the phase modulation of the flat- 
space part, 

2 . O(euj) from derivatives of the phase term, and 

3. 0(ek) from the Christoffel symbols. 


We assume that the light has a much shorter wavelength 
than the gravitational wave (that is, uj >> fc), which is 
surely valid for any real or contemplated detector. We 
further assume that the net phase shift from Eq. ( |40| ) is 
small, meaning euo/k <C 1. Accordingly, we have 

uj euS 2 /k a j ek. (42) 


Hence, to leading order, the electric and magnetic field 
will receive the same phase modulation as the four- 
potential. In fact, all covariant first derivatives of the 
potential receive the same modulation as the potential. 
This also means that our solution respects the Lorentz 
gauge condition. 

In other words, the phase modulation (40) is the 
leading-order effect of the gravitational wave on the elec- 
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tromagnetic field. If we now write 

uj 2 x -uj 2 =uj 2 sm 2 0 cos 20, uj z = uj cos#, z = rcosO 

(43) 

we obtain the expression 0 for the phase. The angular 
dependence we find is confirmed by gravitational wave 
antenna theory (see Maggiore (2008) cf. Eq 9.136 j30l). 


D. Pulses 


The solution (39) corresponds to a single plane wave. 


Just like in Minkowski space, the linearity of the elec¬ 
tromagnetic field equations in a curved space-time (34) 


means we may stack solutions corresponding to what¬ 
ever the light source does. We may replace (39) with 


its Fourier transformation - i.e. an integral over all fre¬ 
quency amplitudes. Consider for example, an experiment 
in which a beam is pulsed at regular intervals. This may 
be modeled as an isolated portion of light, emitted at a 
definite time from the emitter, which travels down the 
cavity to arrive at the observer and, depending on the 
phase, will be observed to arrive either before or after a 
Minkowski photon. If we let the pulse take on the Gaus¬ 
sian wave-packet form at the emitter, in light-frequency 
space it is 


A (w) = exp 


—a (l d — cjo) 


2 


(44) 


and in real space 


A (x M ) ~ cos [cjo (r — t + 4>)] exp 



t + $) 2 , 


(45) 

with <f> given by ( |40| ). This traveling wave packet’s mo¬ 
tion is modulated by the gravitational wave, affecting 
its arrival time at the observer. Because each frequency 
component interacts with the gravitational wave differ¬ 
ently 0 , the pulse shape changes as it travels. 


V. DISCUSSION 


Figure [^summarizes, in mind-map style, the relation¬ 
ship between concepts and their progression as explored 
in this paper. The essential results are Eq. 0 . or Eq. 
(f23l) for the phase modulation, and Eqs. (fTTI) or (|45|) for 


time delays. The well-known schematic picture of a grav¬ 
itational wave detector in terms of a ring of freely falling 
test particles is avoided, and the conundrums associated 
with that picture nn mi do not appear. 

As an alternative picture, we imagine listening to a bird 
song at a rock concert: the bass thump produces local 


changes in the sound speed, inducing a pseudo-Doppler 
shift in the bird song. Or, in our case, the gravitational 
wave modulates space-time’s refractive index. 

Another correspondence is possible with Alcubierre’s 
warp drive [33] space-time, consisting of a time- 
dependent metric for a kind of cavity which skims atop 
a Minkowski background. The respective local contrac¬ 
tion and expansion of space-time in front of and behind 
the moving cavity permit the special relativity-abiding 
passengers of the bubble to experience arbitrarily short 
travel times [34]. The Alcubierre approach to faster-than- 
light travel is not unlike the scenario we have described 
in this paper. The waving space-time provides regions 
of time-dependent contraction and expansion. From the 
wave-optics perspective, the phase velocity (39) modu¬ 
lates about c. A light pulse shot down a waving space- 
time at just the right time can arrive at its destination 
earlier than if the space-time were Minkowski - i.e., if no 
wave were present. An interesting outcome of the geo¬ 
metric optics approach (TT|) is that the expression for the 
time delay holds even when the trajectory is not null, i.e. 
it generalizes trivially to massive bodies. In other words, 
gravitational waves allow for travel times very slightly 
shorter than what is permitted by special relativity, even 
for spaceships. 

The simple form of the solutions we present mean that 
they are suitable for illustrations as animations. In the 
Supplemental Material to this paper we include two user- 
interactive animations written in Python using the mat- 
plotlib library [35]. They animate the time evolution 
of the electromagnetic field as it travels down a cavity 
in a “’waving space-time”. One animation represents 
monochromatic laser light (phase shifts), and the other 
a pulsed beam (time delays). 

Solving the EMW equation under GW modulation (37) 
under more exotic boundary conditions would be intrigu¬ 
ing: for example examining how a waveguide’s propa¬ 
gation modes are altered by the gravitational wave. It 
would be interesting to apply the phase-modulation pic¬ 
ture to extend time-delay interferometry, with unequal 
interferometer arms, moving sources and all [363? to the 
free spectral range. It is possible that phase locking of 
the detector cavity gives rise to resonances, which future 
instruments could exploit. 
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FIG. 3. A mind map showing the relationships between some of the concepts discussed in this paper. There are two routes 
each to the two observables: time delay, and phase modulation. 
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